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Abstract 

Background: Targeted Induced Loci Lesions IN Genomes (TILLING) is increasingly being used to generate and 
identify mutations in target genes of crop genomes. TILLING populations of several thousand lines have been 
generated in a number of crop species including Brassica rapa. Genetic analysis of mutants identified by TILLING 
requires an efficient, high-throughput and cost effective genotyping method to track the mutations through 
numerous generations. High resolution melt (HRM) analysis has been used in a number of systems to identify 
single nucleotide polymorphisms (SNPs) and insertion/deletions (IN/DELs) enabling the genotyping of different 
types of samples. HRM is ideally suited to high-throughput genotyping of multiple TILLING mutants in complex 
crop genomes. To date it has been used to identify mutants and genotype single mutations. The aim of this study 
was to determine if HRM can facilitate downstream analysis of multiple mutant lines identified by TILLING in order 
to characterise allelic series of EMS induced mutations in target genes across a number of generations in complex 
crop genomes. 

Results: We demonstrate that HRM can be used to genotype allelic series of mutations in two genes, BroA.CAXla 
and BroA.METlo in Brossico rapo. We analysed 12 mutations in BroA.CAXlo and five in BroA.METl.a over two 
generations including a back-cross to the wild-type. Using a commercially available HRM kit and the 
Lightscanner™ system we were able to detect mutations in heterozygous and homozygous states for both genes. 

Conclusions: Using HRM genotyping on TILLING derived mutants, it is possible to generate an allelic series of 
mutations within multiple target genes rapidly. Lines suitable for phenotypic analysis can be isolated approximately 
8-9 months (3 generations) from receiving M3 seed of Brossico ropo from the RevGenUK TILLING service. 
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Background 

The identification and analysis of gene mutations in 
plants is fundamental to the investigation of gene func- 
tion. One approach, Targeted Induced Loci Lesions IN 
Genomes (TILLING) was originally developed in Arabi- 
dopsis [1] and has subsequently been successful in a 
range of crop plants [2,3]. In this reverse genetic 
approach, an ethyl methane sulfonate (EMS) mutagen- 
ised population is screened for SNPs within target genes 
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[4]. EMS mutagenesis generates multiple alleles within 
each gene, including nonsense, missense, splicing and 
cis-regulatory mutants, in comparison to T-DNA and 
transposon mutagenesis that generate only knockout 
mutants [5,6]. Analysis of an allelic series can provide 
information on important domains or amino acids 
within the protein of interest. A number of TILLING 
populations have been developed for a variety of crops 
of different genome size, including rice [2], wheat [7,8], 
Brassica rapa [3], Brassica napus [9], Lotus japonicus 
[10], Medicago trunculata [11], Arachis hypogaea [12], 
and Solanum lycopersicum [13]. These populations have 
subsequently been used to isolate mutations in a variety 



o 



© 201 1 Lochlainn et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BIoIVIGCI C6ntr3l commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Lochlainn et al. Plant Methods 201 1, 7:43 
http://www.plantmethods.eom/content/7/1/43 



Page 2 of 9 



of genes, including those involved in starch metabolism 
in Lotus japonicas [14], peanut allergens [12] and a gene 
encoding a fatty acid elongase [9]. 

High resolution melt (HRM) analysis is a technique 
that measures the disassociation of double-stranded 
DNA at high temperature resolution, and permits the 
analysis of genetic variations (SNPs, mutations, methyla- 
tion) in PGR amplicons [15-17]. This technique allows 
genotyping and mutation scanning without the need for 
costly labeled probes, as it uses high fidelity heterodu- 
plex-detecting double-stranded DNA binding dyes, such 
as EvaGreen which exhibit equal binding affinities for 
GG-rich and AT-rich regions and no sequence prefer- 
ence [18,19]. The heteroduplex products are detected by 
the presence of a second low-temperature melting tran- 
sition [20]. This enables the reaction to be performed in 
a single tube, making it cost-effective and suitable for 
high-throughput screening [17]. HRM has been used 
extensively in the genotyping of human tissue samples 
for the identification of genes associated with diseases 
[21,22] as well as identification of clinically important 
fungal species [23]. HRM has been used for quantitative 
detection of adulteration with related species of pants 
used for medicinal purposes [24]. HRM has been used 
previously to identify mutations in TILLING popula- 
tions of tomato [13], wheat [8] and Medaka [25]. In the 
tomato study, HRM was used to identify mutations 
using DNA pools from an EMS mutagenised population 
of plants. The wheat studies demonstrated that HRM 
can be used to isolate novel mutations in starch branch- 
ing enzyme Ila (SBEIIa) genes from an EMS population 
[8] and in mutation scans in mixed PGR amplicons con- 
taining three homeologous gene fragments [19]. It was 
also demonstrated that homozygous and heterozygous 
individuals of a single mutation could be genotyped 
using HRM [19]. 

The genus Brassica includes the closest crop relatives 
of Arabidopsis thaliana, such as B, rapa (A-genome, 2n 
= 2x = 20, -550 Mbp), which includes vegetable crops 
(e.g. turnip, Ghinese cabbage) and oil-seed crops, B. 
oleracea (G-genome, 2n = 2x = 18, > 600 Mbp) which 
includes vegetable crops (cauliflower, broccoli, cabbage) 
and the amphidiploid B, napus (AG-genome, 2n = 4x = 
38, -1100 Mbp), which includes oil-seed crops (canola, 
oilseed rape) and swede. As with many crop plants Bras- 
sica genomes are complex arising from a series of dupli- 
cation events that has resulted in most genes being 
present in multiple paralogous and homeologous copies 
[26,27]. A TILLING population has recently been gener- 
ated in B. rapa to enable gene function to be analysed 
[3] and has previously been used to identify mutations 
in the B.rapa orthologue of the INDEHISGENT gene 
[28]. The aim of this study is to determine if HRM ana- 
lysis can be used to genotype TILLING mutants in the 



crop species B, rapa and thereby to improve the effi- 
ciency of functional genomics approaches in this species. 

Results and Discussion 

Gene selection and TILLING 

In order to demonstrate that HRM is suitable for effi- 
cient genotyping of mutant lines identified using TIL- 
LING, two target genes were chosen, BraA.CAXl.a and 
BraAMETl.a. In A. thaliana, CAXl has been shown to 
encode a calcium-proton antiporter [29] and METl has 
been shown to encode a DNA methyltransferase [30]. 
The CAXl sequence from A, thaliana [Genbank: 
NM_129373.3] was used to identify the homologous 
sequence from B, rapa and from this sequence, 1.5 kb 
including the transcriptional start was used as the target 
for TILLING. To isolate BraAMETl.a from B. rapa R- 
0-18, primers were designed based on the B, rapa 
Ghiifu-401 sequence [Genbank:AB251937] [31]. These 
produced a 2110 bp DNA fragment covering 136 bp 
upstream of the translational start site together with the 
first intron and the first 407 bp of the second exon. 
Both sequences were submitted to the RevGenUK TIL- 
LING service, which resulted in 20 and 17 mutations 
being identified in BraA.CAXl.a and BraA.METl.a 
respectively (Table 1). These mutations included silent, 
missense, nonsense and non-coding mutations. 

Genotyping of TILLING lines using HRM 

From the lines identified using TILLING, twelve 
[BraA.caxl.a-l, -4, -6, -7, -8, -10, -11, -12, -14, -16, 
-18, -20) and five {BraA.metl.a-3, -6, -7, -10 and -14) 
lines carrying mutations in BraA.CAXl and BraA. 
METl. a respectively were chosen to be characterised 
based on the nature and position of the mutation. In 
order to genotype the lines, a HRM assay was devel- 
oped for each gene. To increase the specificity of the 
HRM reaction a nested PGR approach was used for 
the BraA.caxl.a lines, first using the TILLING primers 
to amplify genomic DNA, prior to using the HRM pri- 
mers with this product. In total, three primer pairs 
were designed to genotype all 12 mutations within the 
BraA.CAXl.a gene (amplifying fragments of 200, 203 
and 260 bp; Table 2) and one primer pair was designed 
to genotype all the mutations within the BraA.METl.a 
gene (amplifying a 399 bp fragment; Table 2). Indivi- 
dual M3 plants of the TILLING lines and wild-type R- 
0-I8 plants were then genotyped using the HRM assay 
using a single primer pair per plate and three technical 
replicates per plant, allowing multiple alleles to be gen- 
otyped in a single scanning run. Analysis of the differ- 
ence plots for the fluorescent signals detected 
individual plants that were wild-type, homozygous or 
heterozygous for the mutation in alleles of both genes 
(Figure 1). Three distinct groups can readily be 
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Table 1 Nucleotide changes in BraA.cax1.a and BraA.metl.a TILLING mutants. 



Line^ Mutation Change'' Position (base-pair)*^ Position (amino-acid)^ Type of mutation^ 



BraA.coxl .a-1 


G to A 


881 




Non-coding 


BroA.coxl .0-2 


C to T 


860 




Non-coding 


BroA.coxl .0-3 


G to A 


545 


2 


Non-coding 


BroA.coxl .0-4 


G to A 


771 


77 


Missense 


BroA.coxl .0-5 


C to T 


494 




Non-coding 


BroA.coxl .0-6 


G to A 


517 




Non-coding 


BroA.coxl .0-7 


G to A 


673 


44 


Missense 


BroA.coxl .0-8 


G to A 


599 


1 9 


Missense 


BroA.coxl .0-9 


G to A 


645 


35 


Missense 


BroA.coxl .0-1 0 


C to T 


736 


65 


Missense 


BroA.coxl .0-1 1 


G to A 


795 


85 


Missense 


BroA.coxl .0-12 


C to T 


708 


56 


Missense 


BroA.coxl .0-1 3 


C to T 


689 


49 


Silent 


BroA.coxl .0-14 


C to T 


661 


40 


Missense 


BroA.coxl .0-1 5 


C to T 


658 


39 


Missense 


BroA.coxl .0-1 6 


C to T 


421 




Non-coding 


BroA.coxl .0-1 7 


G to A 


926 




Non-coding 


BroA.coxl .0-1 8 


C to T 


81 8 


92 


Silent 


BroA.coxl .0-19 


C to T 


732 


64 


Missense 


BroA.coxl .0-20 


C to T 


690 


50 


Nonsense 


BroA.metl .0-1 


C to T 


3516 


248 


Silent 


BroA.metl .0-2 


G to A 


3334 


1 88 


Missense 


BroA.metl .0-3 


G to A 


3334 


1 88 


Missense 


BroA.metl .0-4 


G to A 


3084 


1 04 


Silent 


BroA.metl .0-5 


C to T 


3081 


1 03 


Silent 


BroA.metl .0-6 


G to A 


3252 


243 


Nonsense 


BroA.metl .0-7 


G to A 


3500 


243 


Missense 


BroA.metl .0-8 


G to A 


301 2 


80 


Silent 


RrnA mptl n-Q 
Dl Ur\.i 1 IcL 1 .U y 


G to A 




241 


^ilpnt 
jllcl 1 L 


BroA.metl. 0-1 0 


CtoT 


3335 


188 


Missense 


BroA.metl. 0-1 1 


CtoT 


3069 


99 


Silent 


BroA.metl. 0-1 2 


G to A 


3330 


186 


Silent 


BroA.metl. 0-1 3 


CtoT 


3324 


184 


Silent 


BroA.metl. 0-1 4 


G to A 


3213 


147 


Nonsense 


BroA.metl. 0-1 5 


CtoT 


3265 


165 


Silent 


BroA.metl. 0-1 6 


G to A 


3330 


186 


Silent 


BroA.metl. 0-1 7 


G to A 


3258 


162 


Silent 



^ - TILLING line, ^ - Nucleotide change in mutant line, ^ - Position of nucleotide change in base-pairs, ^ - Position of amino-acid change in peptide, ^ - Type of 
mutation caused by nucleotide change. Non-coding - nucleotide change not within coding region, missense - nucleotide change results in an amino-acid 
change, nonsense - nucleotide change results in the introduction of a premature stop codon, silent - no change in amino-acid. 



observed in the difference plots. To confirm these 
results and assign zygosity groups, the PGR products 
from selected lines were sequenced. Sequencing the 
HRM PGR products identified by HRM as heterozy- 
gous from BraA.caxl.a-11 and BraA.metl.a-6, showed 
a double peak at the position of the mutation (Figure 
2), confirming that these were indeed heterozygous for 



the mutations. In addition, sequencing of PGR pro- 
ducts from BraA.caxl.a-1 and BraA.metl.a-6 identified 
as homozygous by HRM, confirmed that they were 
homozygous for the mutation (Figure 2). These results 
demonstrate that HRM can be used both for rapid 
genotyping of TILLING mutants and for distinguishing 
wild-type, heterozygous and homozygous plants. The 
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Table 2 Details of primers used for gene Isolation and HRM analysis. 



Primer^ 


Gene'' 


Stage*^ 


Mutation Targeted^ 


Sequence 5' to 3'^ 


DrCAXl MLUNG_l-orwara 


braA.LAXl .a 


TILLING 




hk3l\Lyl\ 1 II CC 1 ACCCA 1 Ci 1 C 


DrCAXl MLLING_Keverse 


A r~ A\/ 'i ^ 

braA.LAXl .a 


TILLING 




r~ r~ A r~ r~ r~r — r k k \ \ r i i i i a Tt — vr~ 

CCACCCC 1 AA 1 1 G 1 1 1 1 A 1 G 1 G 


BrCAXl HKI\/ll_l-orwara 


timA.LAXi .a 


HRM 


6,1 6 


Tr~f — vr~r~ a a t — rT/^r~f — rr — ri^ at 

1 CC 1 CGAAG 1 1 GCC 1 C 1 GA 1 


BrCAXl H KM l_Ke verse 


broA.LAXl .0 


HRM 


6,1 6 


GC 1 GC 1 GACCA 1 1 G 1 1 CC 1 G 


BrCAXl HKI\/li_l-orwara 


n y^ A r~ A\/ 1 ^ 

braA.LAXl .a 


HRM 


7,8,1 2,14,20 


LjLALjLjAALAA 1 GG 1 CAGCAG 


BrCAXl HKI\/li_Keverse 


n y^ A r~ A\/ 1 ^ 

braA.LAXl .a 


HRM 


7,8,1 2,14,20 


CGAGAA 1 GAG II C II GGAGA 1 1 


BrCAXl HKI\/lD_l-orwara 


n y^ A r~ A\/ 1 ^ 

braA.LAXl .a 


HRM 


4,10,1 1,18 


\ CCAGAAGG 1 1 CCA 1 ACAAAGG 


DrCAXl HKlVlD_Ke verse 


braA.LAXl .a 


HRM 


/I 1 r\ 1 1 1 o 

4,1 0,1 1 ,1 8 


r~ c ( — V( — m — V( — 1 1 1 A c — r a a T( — r/ — vrc 
CGC 1 G 1 1 C 1 G 1 1 1 AG 1 AA 1 G 1 G 1 1 G 


Di IVlcL 1 dr^ 


RrnA MFT1 n 

Dl Ur\.iv\C I 1 .U 


Ocl Ic IbUldLIUI 1 






BrMetlaFS 


BraA.METl.a 


Gene isolation 




CTOACAATGGTGCTGTOTOAG 


BrMetlaF6 


BraA.METl.a 


Gene isolation 




GTCTGGmGATGGCAGAGGACG 


BrMetlaR4 


BraA.METl.a 


Gene isolation 




CGCACTCGCAGCTTCAGAAAC 


BrMetlaRS 


BraA.METl.a 


Gene isolation 




AAATAGAGATOTOAACACCCCC 


BrMetlaR6 


BraA.METl.a 


Gene isolation 




TOGATACCACAGGTGCCATGC 


BrMetlaHRMF 


BraA.METl.a 


HRM 


3,6,7,10,14 


TGCGATGATAAGGAGAAAGG 


BrMetlaHRMR 


BraA.METl.a 


HRM 


3,6,7,10,14 


mGCCGTCTCATCCAAAC 



^- Primer name, "^-gene primer designed to, ^-stage of analysis primer used, ^-line primers used to genotype, ^- sequence of primer 



results also demonstrate that multiple alleles can be 
genotyped using a single primer pair on one PGR 
plate, allowing multiple lines and individual plants to 
be genotyped simultaneously. 



Backcross analysis 

To reduce the mutation load, reciprocal crosses of the 
TILLING lines and wild-type R-o-18 plants were per- 
formed with both homozygous and heterozygous mutant 





82 83 



85 8B 87 88 89 

Temperature °C 
Difference curves 




Temperature °C Temperature °C 

Figure 1 HRM analysis for Brassica rapa R-o-18 wild-type and individual plants from M3 BraA.caxha-ll and BraA.metha-6 TILLING 
mutants. A: Normalised melt curves for fluorescent signals from DNA strand dissociation of triplicated technical replicates of B. rapa wild-type 
(light green lines) and individual plants of M3 TILLING mutant BraA.caxl .a-1 1 , plant BraA.caxl .a-1 1 A (red), plant BraA.caxl.a-1 IB (red), plant BraA. 
caxl.a-llC (blue), plant BraA.caxl.a-1 ID (light green), plant BraA.caxl.a-1 IE (dark green). B: Normalised melt curves for fluorescent signals from 
DNA strand dissociation of triplicated technical replicates of B. rapa wild type (red lines) and individual plants M3 TILLING mutant BraA.metl .a-6, 
plant BraA.metl .a-6-1 (light blue), plant BraA.metl. a-6-2 (green), plant BraA.metl. a-6-3 (orange), plant BraA.metl .a-6-4 (purple), plant BraA.metl. a- 
14-1 (blue). C: Difference plot of genotypes' fluorescence normalised to wild type samples (light green lines, i) and individual plants of M3 
TILLING mutant BraA.caxl.a-1 1 , plant BraA.caxl .a-1 1 A (red, ii), plant BraA.caxl.a-1 IB (red, ill), plant BraA.caxl.a-1 IC (blue, iv), plant BraA.caxl.a-1 ID 
(light green, v), plant BraA.caxl.a-1 IE (dark green, vi). D: Difference plot of genotypes' fluorescence normalised to wild type (red, i) and individual 
plants of the M3 TILLING mutant BraA.metl. a-2, plant BraA.metl. a-6-1 (light blue, ii), plant BraA.metl. a-6-2 (green, ill), plant BraA.metl. a-6-3 
(orange, iv), plant BraA.metl. a-6-4 (purple, v), plant BraA.metl. a-1 4-1 (blue, vi). 
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Ac * T * C C T 0 C C 










1 1 ~ — 

I'I'IMMMa'I 


7^ ^ 

Qc GTG GTTT G 


— 

TTCTC G*CG G 


Figure 2 Sequencing chromatograms of M3 generation of B. 
rapa R-o-18 mutants. A: BraA.coxl .a-ll , gene sequence from 788 
bp to 807 bp shown. B: BraA.metl.o-6, gene sequence from 3246 bp 
to 3265 bp shown. C: BroA.coxl.o-1, gene sequence from 868 bp to 
888 bp shown. D: BraA.metl.o-6, gene sequence from 3246 bp to 
3265 bp shown. Rectangular box indicates zygosity state at the site 
of EMS induced transitional mutation. 



plants. A heterozygous plant from the BraA.metl.a-6 
line was crossed to a wild-type plant and ten Fi plants 
were genotyped. The results from the genotyping 
showed the expected segregation of the mutant allele, 
with 5 plants having a heterozygous mutation and 5 
plants being wild-type for the mutation. The genotyping 
results for the plants from crosses of the BraAxaxl.a 
lines also showed the expected genotyping results (Table 
3). For example, the Fi plants of a cross using a homo- 
zygous BraA.caxl.a-1 plant were all heterozygous for 
the mutation (n = 5). In contrast, the Fi progeny of het- 
erozygous plants showed segregation for the mutation 
(Table 3). These results demonstrate that HRM can be 
used to follow mutations and zygosity through multiple 
generations and crossing events, making it an ideal 
method for genotyping of TILLING mutants. 



Conclusions 

To analyse an allelic series of mutants generated using 
the TILLING approach efficiently, a cost effective, high- 
throughput genotyping technique is required. Techni- 
ques that have been used previously include PCR-RFLP, 
but this relies on the mutation causing a change in a 
restriction site [32] and TaqMan assays, which can 
require expensive probes and assay optimisation [18,33]. 
Alternatively, PGR products covering the mutations 
could be sequenced. However, this is time consuming 
and potentially more expensive; HRM analysis of sam- 
ples in triplicate costs -30% of Sanger sequencing costs. 
The results presented here demonstrate that HRM can 
be used to genotype mutations identified by TILLING 
from EMS populations in B. rapa. The results also 
demonstrate that a single PGR primer pair could be 
used to identify mutations at multiple positions with the 
PGR fragment. This reduces the number of primer pairs 
needed to genotype an allelic series and hence the cost. 
Due to the nature of the HRM process, the DNA extrac- 
tion, PGRs and melt curve analysis can all be performed 
in 96-well plates. This enables high-throughput genotyp- 
ing of multiple alleles of a number of genes to be per- 
formed quickly and efficiently. This is particularly 
important in complex crop species such as B rapa, 
where paralogous genes may be present and accumula- 
tion of mutations may be required to generate a 
phenotype. 

Methods 

Plant material 

Seed of the highly inbred, homozygous, self-compatible, 
rapid cycling B, rapa ssp. trilocularis line R-o-18 was 
used as wild-type control. Seeds for lines (M3 genera- 
tion) containing mutations in target genes were obtained 
from the R-o-18 TILLING population through the 
RevGenUK service [3,34]. 

Plant growth conditions 

Plants were grown under glasshouse conditions with a 
16 h photoperiod at 22.3 ± 4°G and 13.3 ± 2°G mean 
day and night temperatures respectively. Seeds were 
sown directly into pots (height 11.3 cm; diameter 13 
cm) containing Levington M3 high nutrient compost 
(Monro Group, Ghichester, UK) and irrigated twice 
daily. Following three weeks of vegetative growth. Vita- 
feed® 2-1-4 nutrient solutions (N-P-K: 16-8-32 + micro- 
nutrients; Vitax Ltd., Goalville, Leicestershire, UK) were 
applied to plants weekly at a rate of 3 g 1'^ via overhead 
irrigation. Plants initiated flowering after approximately 
6 weeks of growth and seeds were collected approxi- 
mately 4 weeks thereafter. Perforated bread bags (150 
mm X 700 mm, WR Wright & Sons Ltd, Liverpool, UK) 
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Table 3 HRM genotyping results of progeny of BraA.cax1,a and BraA.metl.a lines crossed to wild-type. 



Line^ 


Plant"^ 


Number of heterozygous plants*^ 


Number of wild-type plants'' 


BraA.CQxl .a-1 


E 


5 


0 


BroA.coxl .0-4 


C 


2 


0 


BroA.coxl .0-6 


B 


4 


2 


BroA.coxl .0-7 


A 


3 


3 


Dl Ur\.LUA 1 .U O 


A 


A 


1 
1 


BroA.coxl.o-10 


A 


3 


2 


BroA.coxl. 0-11 


B 


3 


3 


BroA.coxl. 0-1 2 


A 


3 


3 


BroA.coxl. 0-1 4 


B 


3 


0 


BroA.coxl. 0-1 8 


D 


2 


4 


BroA.coxl .0-20 


A 


2 


2 


BroA.metl.o-6 


A 


5 


5 



^ - TILLING line crossed to B. rapa R-o-18 wild-type, ^ - Individual IVI3 plant identifier' ^ - Number of individual Fi heterozygous plants identified by HRM' ^ - 



Number of individual Fi wild-type plants identified by HRM 



were used to enclose inflorescences to prevent cross- 
pollination from neighboring plants. 

In silico data mining to identify gene sequences 

Basic Local Alignment Search Tool (BLAST) [35] soft- 
ware was employed to identify sequences from all Bras- 
sica spp. that were orthologous to the A. thaliana CAXl 
gene [Genbank: NM_129373.3]. AU A. thaliana 
sequences of interest were retrieved from GenBank [36]. 
For the BLASTn algorithm [37] the entire GenBank 
database (March 2009) was interrogated to identify all 
Brassica sequences sharing > 80% sequence identity at 
an e-value of < le'^^ with the A. thaliana query 
sequence. The Brassica Genome Sequencing Project 
(^rGSP) database [38] was interrogated (March 2009) 
using the WuBLAST [39] algorithm for all Brassica 
BAC-end sequences, fully sequenced BACs, ESTs and 
genome survey sequence which shared > 80% sequence 
identity with this query sequence. A. thaliana genes 
sharing high sequence similarities to the query gene Le, 
members of the same gene family including AtCAX2 
[Genbank:NM_001036119.1] and A/:CAX3 [Genbank: 
NM_1 15045.3] were used to interrogate both databases. 
After in silico comparative sequence analyses, Brassica 
sequences were assigned to the member of the A. thali- 
ana gene family with which they shared the highest 
sequence identity, thus removing false positives. All 
selected sequences which were orthologous to AtCAXl 
were formatted into FASTA and entered into a contig 
assembly program (ContigExpress, VectorNTI 11, Invi- 
trogen. Paisley, UK) using default settings to form 
longer contiguous Brassica sequences and identify possi- 
ble locus specific paralogues. These nascent Brassica 
gene sequences were aligned to the A, thaliana gene of 
interest, flanked at either end with 1 Mbp of sequence. 



using AlignX software at default settings. The genomic 
structure of these contigs was elucidated and annotated 
in the Vector map software (VectorNTI 11). Primers 
(Table 2) were designed to amplify ~ 1.2 kbp fragments 
of BraA.CAXl.a beginning upstream from the deduced 
transcriptional start site to a region within the gene 
using Primer 3 (Version 0.4.0) [40] with the parameters 
"Max Self Complementarity" and "Max 3' Self Comple- 
mentarity" adjusted to 2, to avoid hairpin loops and 
potential dimerisation. Primers to amplify the BraA. 
METl.a gene (Table 2) from B. rapa R-o-18 were 
designed based on the sequence described by Fujimoto 
et al. [[31], Genbank:AB251937]. 

Amplification of locus specific sequences from B. rapa R- 
0-18 

DNA was extracted from leaves of wild-type B. rapa R- 
0-I8 using the GenElute® Plant Genomic DNA Mini- 
prep kit (Sigma- Aldrich, GilUngham, UK) according to 
the manufacturer's instructions and subjected to PCR 
amplification using primers (Table 2) designed to Bras- 
sica contig sequence. Amplified fragments were ligated 
in plasmid pCR8®/GW/TOPO® (Invitrogen, Paisley, 
UK) and subsequently sequenced. 

Plant genomic DNA extraction 

Genomic DNA was extracted from the BraA.metl.a 
plants using the DNeasy Plant Min kit (Qiagen, Crawley, 
UK) following the manufacturer s instructions or using a 
rapid salt-extraction method [41]. For the BraA.caxl.a 
plants, a high throughput extraction method was devel- 
oped based on a method developed by Dr. M. Poole 
(pers. comm.). A 0.25 cm^ diameter leaf disc from a 
developing or mid-expansion leaf from each plant was 
placed in a 96-well plate (Qiagen Collection Microtubes, 
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Qiagen, Crawley, UK), with a single 3 mm stainless steel 
ball and a 600 [d aliquot of DNA extraction buffer (200 
mM Tris-HCL pH 7.5, 250 mM NaCl, 25 mM EDTA 
and 0.5% w/v SDS) added to each well. Suspensions 
were shaken for 3 min in a ball mill (8000 Mixer Mill, 
Glen Creston, Stanmore, UK) until the majority of tissue 
was disrupted. Plates were subsequently centrifuged in a 
plate centrifuge (Centrifuge 5804, Eppendorf UK, Cam- 
bridge, UK) at 5,600 g for 20 min. After centrifugation 
300 \i\ of the supernatant was transferred to a fresh 
plate (Thermofisher storage plate with rubber sealing 
mat. Fisher, Loughborough, UK), containing 300 [i\ of 
100% isopropanol, mixed and incubated for 3 min at 
room temperature. Plates were then centrifuged at 5600 
g for 20 min before discarding the supernatant. A 100 
[A aliquot of 70% (v/v) ethanol was added and mixed by 
shaking vigorously for 5 min, centrifuged at 5,600 g for 
20 min. The supernatant was discarded and dried in a 
vacuum manifold for 10 min. DNA was resuspended in 
20 \i\ sterile H2O. The DNA concentration was not 
quantified and standardized due to the number of sam- 
ples extracted simultaneously. 

High Resolution Melt (HRM) analysis 

Primers were designed for HRM analysis to amplify 100- 
250 bp regions within the BrciA..CA.Xl ,ci gene and a 399 
bp region of BraAMETl.a, which were originally used 
to identify mutations by TILLING. Three sets of HRM 
primer pairs were designed for BraA,CAXl,a and one 
for BraA.METl.a (Table 2), using Primer 3 (Version 
0.4.0) with the parameters "Max Self Complementarity" 
and "Max 3' Self Complementarity" adjusted to 2, and a 
Tm > 60°C, to avoid hairpin loops and dimerisations. 
Primer sequences were examined further by primer 
thermodynamic software (Vector NTI 11) to confirm 
that no secondary structures were likely to form. HRM 
primer amplification efficiencies and specificities, were 
determined in vitro by PCR amplification of wild-type B. 
rapa R-o-18 DNA in reaction volumes containing 10 \A 
of GoTaq buffer (5x), 5 \A of MgCl2 (25 mM), 1 \A of 
dNTPs (10 mM), 5 (il of both forward and reverse pri- 
mers (10 pmol (il'^), DNA (10 ng), 0.2 \A of Taq poly- 
merase (GoTaq) and sterile de-ionised water (SDW) to 
a total volume of 50 \A, PCR conditions were 95°C for 5 
min and 30 cycles of 95°C for 10 sec, 55°C for 10 sec 
and 72°C for 30 sec. 

Nested PCR approaches to amplify locus specific 
fragments for HRM analysis 

The nested PCR approach was used for the BraACAXl, 
a lines. DNA extracted from TILLING lines were ampli- 
fied with TILLING primers using Phusion™ high fide- 
lity DNA polymerase (New England Biolabs, Hitchin, 
UK) to ensure sequence integrity of amplified sequences. 



Reaction volumes comprised of 4 \A of Phusion HF (5x), 
0.4 (il of dNTPs (10 mM), 1.4 \A of forward and reverse 
BraA.CAXl.a TILLING primers (10 pmol (iL^), 1 (il 
DNA, 0.6 \A of DMSO, 0.2 \A of Phusion DNA polymer- 
ase and SDW to a total volume of 20 (il. PCR conditions 
involved an initial denaturation at 98°C (30 sec) followed 
by 30 cycles of denaturation at 98°C (5-10 sec), anneal- 
ing at 58°C (15 sec) and extension at 72°C (30 sec). 

HRM reactions 

HRM were performed using Type-it® HRM™ PCR kit 
(Qiagen, Crawley, UK) following manufacturer's instruc- 
tions; 1 (il of the initial nested PCR reaction for the 
BraA.caxl.a lines and 1 \\\ of genomic DNA for the 
BraA.metl.a lines was added as template to the HRM 
PCR reaction mix comprising 5 (il of 2 x HRM PCR 
Master Mix (including EvaGreen fluorescent dye), 1.75 
(il of forward and reverse HRM primers (10 pmol (iL^) 
and SDW to 10 \A, Using an Eppendorf 96-well Master- 
Cycler 5331, reaction volumes were subjected to 95°C (5 
min), followed by 30 cycles of 95°C (10 sec), 55°C 
annealing (10 sec) and 72°C extension (30 sec). Included 
within each 96-well analysis plate were samples in tripli- 
cate, also triplicate wild-type DNA samples and two sets 
of negative controls, one without DNA added and one 
without HRM primer added. 

Lightscanner^M analysis of HRM products 

PCR plates containing HRM amplicons were placed into 
a Lightscanner™ (Idaho Technology, Salt Lake City, 
USA) and the temperature raised to 60°C for 5 mins to 
ensure all samples were equilibrated. The melt tempera- 
ture range was set at 63.3°C-95.4°C with a ramp setting 
at 0.1°C and a second hold at each step. Exposure was 
set at "Auto", background correction to exponential, 
curve shift to 0.020, standards to 'Auto Group', and sen- 
sitivity at normal +2.8. HRM analysis was then per- 
formed on the dissociation of double-stranded DNA 
PCR products, which had been saturated with the low 
PCR-toxic dye, EvaGreen from the initial HRM PCR. 
The LightScanner™ Data Analysis software (Version 
2.0, Idaho Technology) was used to analyse the data and 
produce normalised disassociation curves and difference 
plots. 
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